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I. INTRODUCTION 


The analysis presented in this report is concerned with 
the effect of residual stresses on the life-lengthening of 
structural members. In particular, it has been observed that 
in specimens subjected to compressive residual stresses, which 
occur as a result of plastic deformation at a stress concen-~ 
tration, the life of the specimen is extended (Ref. 1-4). 
However, in some later tests involving both load and. time 
effects, the residual compressive stresses have been shown 
to shorten the specimen lifetime (Ref. 5). This apparent 
contradiction is no doubt due to the complex phenomenological 
behavior of the material which is not very well understood 
at the present time. | 

To better understand the material behavior as discussed 
briefly above (and referred to as load-time dependent relaxa- 
tion of residual stresses) an analytical approach has been 
developed to complement the existing experimental approaches. 
The techniques involved and the results of this analysis are 
presented in Section II. The conclusions and recommendations 
are presented in Section III. The supporting documentation 


for this analysis is included in the Appendices. 
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DISCUSSION OF RESULTS 

The test specimen made of aluminum 7075 T651 (bare) and 
shown in Figure 1 has been used in this analysis. Creep test 
data for this specimen geometry and material have been made 
available to ASIAC (Ref. 6). To analyze this specimen the 
computer code CREEP-PLAST (Ref. 7) has been used. CREEP-PLAST 
is a two-dimensional finite element code for the creep and 
plasticity analysis of metal structures operating at high 
temperatures. | | 

CREEP-PLAST considers instantaneous time-independent 
elastic-plastic and time-dependent creep simultaneously in 
an incremental procedure which requires the user to select 
his own discrete time and/or load increments. The method of 
calculation is based on the assumption that the total strain 
at any instant consists of elastic, plastic, and creep parts 
with essentially no interdependence between the plastic and 
creep components. The program is capable of predicting de- 
formations, stresses, and total accumulated strains as func- 
tions of time and/or load for user specified thermal and 
mechanical loadings which can be varied independently. 

A more detailed discussion of the capabilities of the 
CREEP-PLAST code can be found in Ref. 7. 

For the prediction of creep strains, the user must first 
specify the creep strain as a function of the three parameters_ 


Ex+>s €, and r, which are then used in the following relation: 


ea. = e,(l-e**) + et 
where 
es. limiting primary creep strain 
r = exponential term in primary creep enn 
a = secondary creep strain rate 
E, = creep strain 
t = time 
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Once the quantities Eas es and r have been determined for 
each creep test at a different stress level it 1s necessary to 
express these quantities (€,, E and r) as functions of the 
stress level. The actual creep test data reduction and pro- 
cedures required to determine Eas es and r as a function of 
the stress level are explained fully in Appendix A. The imple- 
mentation of these relations just discussed requires that the 
user modify Subroutine SHIFT in CREEP-PLAST. The modifications 
to Subroutine SHIFT are included in Appendix B. 

Also it is necessary to express the tangent modulus as a 
function of the stress level. To do this a representative 
stress-strain curve from MIL-HDBK-5 has been used. This 


stress-strain curve is then fitted with a polynomial of the 


form 


eS. vas sk DO + co? + age 


from which the tangent modulus CE.) is derived. 


a E = b + 2co0 + 3do 

The actual procedure is discussed more fully in Appendix B. 
Also in Appendix B is a listing of Subroutine MATCON which 
the user must modify by incorporating the tangent modulus as 
a function of the stress. 

The finite element model (Figure 2) considers only one- 
quarter of the specimen geometry (See Figure 1) since symmetry 
can be used to reduce the problem size. The boundary condi- 
tions, coordinate axis, and loading are also indicated in 
Figure 2. Appendix C contains several enlarged areas of 
Figure 2 along with element and node numbers for the conven- 
ience of the reader. 

The input data echo from CREEP-PLAST is contained in 
Appendix D. Appendix D, along with Appendix C provides a 


complete description of the finite element model. 
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The particular load-time history used in this analysis 
is shown in Figure 3. The temperature was assumed to be a 
constant 85°F throughout the analyses. 

Since the CREEP-PLAST finite element code can handle 
only mono-modulus (that is, either tension or compression, 
but not both) type properties, it was necessary to choose 
either tension or compression creep properties for the analysis. 
The material constants used in the elastic constituitive rela- 
tions (i.e., E, v) are very nearly the same for tension and 
compression, however, as can be seen from Figures A-1 and A-2 
of Appendix A, the creep behavior in tension or compression 
1s quite different. Since the important material behavior 
of load-time dependent relaxation of residual stresses occurs 
as a result of an applied tensile load, it was somewhat 
arbitrarily decided to use the creep properties corresponding 
to the 85°F tension creep tests. In examining the impact of 
this decision on the results of this analysis, it can generally 
be concluded that the significant effect of load-time depen- 
dent relaxation of residual stresses Pel still be observed, 
however, the magnitude of total strain relaxation may be 
different. In particular note that load-time dependent relax- 
ation of residual stresses will occur ina relatively short 
time period and thus the secondary creep term (e£) will be 
negligible. Also note that the primary creep rate (r) is the 
same for tension and compression. Thus the difference between 
tension and compression creep behavior for small time is due 
to the quantity Ey. Thais duentaty elects the amount of 
residual strain relaxation, not the rate. » Thus even though 
the choice of tension creep properties is somewhat arbitrary 
the basic conclusions will remain the same. 

An examination of the first iteration, which is an elas- 
tic solution, shows that the stress concentration factor is 
approximately 2.5, which compares favorably with theoretical 


derivations for the stress concentration factor. 
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A plot of the total strain (TOTAL = ELASTIC + PLASTIC 
+ CREEP) is shown in Figure 4 for t < 87,000 sec. It is 
observed that from about t = 30 sec. to t = 6000 sec., that 
load-time dependent relaxation of residual stress occurs. 
In particular at t = 30 sec. (start of continuously applied 


compressive load) the total (effective) strain is 8310 uw in/in. 


At t = 6000 sec., the effective strain is 2682 yp in/in. From 

t = 6000 sec., the behavior is quite clearly a result of creep 
strain under a compressive load with both primary and secondary 
creep in evidence. 

The CREEP-PLAST code does not print out accumulated creep 
strain, however, shown in Figure 5 is a plot of the elastic 
strain and the non-elastic (plastic and creep) strain for the 
time period involved. After the initial relaxation period 
the accumulated strain is primarily due to creep. Appendix E 


contains the CREEP-PLAST code computer output for selected 


times. 
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CONCLUSIONS AND RECOMMENDATIONS 

The analytical technique and results presented herein 
for the load-time dependent relaxation of residual stresses 
in an aluminum test specimen has been shown to be a valuable 
addition to the experimental work previously done in this 
area. This analysis is essentially a demonstration of the 
value of using an analytical approach in order to gain a 
better understanding of complex material behavior. However, 
as in all efforts of this nature, it is possible to improve 
the results. 


Specifically, the following recommendations are made: 


(1) Determine by test within the same specimen lots the 
actual nonlinear stress-strain curve. Thus, creep test 
data will have an associated stress-strain curve. In 
this manner differences between handbook data (1.e., 
MIL-HDBK-5) and the actual specimen behavior are excluded 


from the analysis. 


(2) Utilize linear regression analysis to minimize any 


curve-fitting errors. 


(3) Consider code modifications to CREEP-PLAST in order to 
handle tension and compression properties. At the moment 
there is not a unified theory that considers the combined 
effects of multiaxial stresses and the effect on creep 
properties, so any modifications at this point would be 


somewhat empirical in nature. 


(4) Consider the possibility of load-time-temperature 
dependent relaxation of residual stresses. Such a 
consideration would involve a larger testing matrix to 


establish property trends as a function of temperature. 
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APPENDIX A 


REDUCTION OF CREEP TEST DATA 
FROM 
ALUMINUM 7075 T651 TEST SPECIMEN 


The creep test data (Ref. 6) used to establish the 
creep strain relations for CREEP-PLAST are shown in Tables A-l 
through A-15. This data is for aluminum (7075 T6591 bare) at 
two different temperatures (85°F and 135°F) and for varying 
stress levels in tension and compression. Since the current 
version of CREEP-PLAST can handle only mono-modulus type 
properties, only the 85°F tension creep test data is reduced 
and utilized in this analysis. 

CREEP-PLAST requires that the creep strains be expressible 


in the following form: 


ba e,(1-e77") + et 
where 
E, = creep strain 
ey = limiting creep strain for primary creep 
r = exponential rate term 
= = secondary creep strain rate 
t = time 


The quantities Ey > E_, and r must be expressed as a func- 
tion of the stress level for use in Subroutines MATCON and 
SHIFT (Sée Appendix B). 

The first step in the data reduction is to plot the 
data of Tables A-1 through A-15. This plotting has been 
shown in Figures A-1 and A-2 along with the fitted curve 
used in CREEP-PLAST. To facilitate the determination of 


the quantities Ei» € and vr, the secondary creep strain rate 


m 
has been estimated visually from Figures A-1 and A-2. By 
choosing two points in the primary creep region--say Ce, ty) 


and Cee »t.)--then two simultaneous equations can be 
2 


used to determine €, and r, thus, 


+ 
_ -“rt ‘ 
E = E. (l-e 1) + ent1 
1 
: -rt : 
os = Ee, (l-e 2) + Et 


from which, by trial and error, r can be found from 


Knowing r, then €, can be found quite easily. 


" 

Although other methods for determining 4s Ee and r can 
be used (i.e., regression analysis to minimize error) this 
approach is sufficiently accurate for the purposes of this 
analysis. 

The quantities E+ > E and r, determined by the above 
approach, are shown in Table A-16. 

Next, it is necessary to express the quantities 


Sus En and vr as functions of stress. This was accomplished 
by assuming that either a straight line or a functional of 
the form 


. 2 3 


€ er = a + bo + co + do 


t? 


can be used to give an adequate fit to the data. The curves, 
along with the data points from Table A-16 are shown in 
Figure A-3. Since the trend for r is not well defined, a 
value of r = 2.0 x 10° */sec was used. The final equations 
are listed below (and are also included in Appendix B) for 
use in subroutines MATCON and SHIFT of the CREEP-PLAST code: 


a 


o < 43,000 psi. E = (6.325 x 10° “)o 
43,000 < o < 66,000 psi. En = ((6.025 x 107”) 
+ (6.569 x 10°7%)(o - 43,000) + (2.668 x 10729) 
7 2 -4 
x Co - 43,000)")(o - 43,000) + 2.591 x 10 
66,000 < o < 71,000 psi. E, = ((3.628 x 0°) 


400903. 10°) (e166 4000) oo L000 a Fubes Gee 


o> 71,000 é& = (1.8 x 10°°)(o - 71,000) + 1.5 x 107 


o < 66,000 €, = (2.045 x 


66,000 < o < 72,000 Ee, 


x (o - 66,000) + (2.041 


uy 


3 


107°)o 


7 


CO208S e110 Ps & C3875 107" 5 


x 1077)(o = 66,000)*)(o - 66,000) + 135 


o >72,000 e, = (2.271 x 107+)(o - 72,000) + 602.1 


t 


for all o> 0 vr = 2.0 x 10 ‘/sec 


TABLE A-1l CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


fe 385 PR STRESS = 30,000 psi. Compressive 
Time (sec) Strain (y in/in) Time (sec) Strain (yp in/in) 

10 4 54000 . 60 
30 8 57600 68 
60 8 | 61200 72 

120 12 64800 76. 
180 12 684.00 80 
240 12 72000 80 
300 12 75600 80 
360 16 79200 8 8 
420 Lb: 82800 88 | 
480 16 86400 88 
540 16 90000 92 
600 5 we 93600 J2 
900 16 97200 J6 
1200 | 16 100800 96. 

1500 20 104400 — 

1800 24 108000 100 
2400 24 111600 104 
3000 ot 115200 104 
3600 28 118800 104 
5400 28 122400 104 
7200 | 28 | 126000 104 
10800 . ~ 39 129600 108 
“14400 ran | 133200 108 
18000 48 136800 108 
21600 48 - 140400 108 
25200 52 144000 112 
28800 52 147600 116 
32400 52 151200 124 
36000 52 154800 128 
39600. | 56 158400 128 
43200 96 165600 © 132 
46800 56 169200 132 


50400 96 | 172800 — OZ 


TABLE A-2 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 
T = 85°F, STRESS = 43,000 psi. Compressive 
Time (sec) Strain (yp in/in) Time (sec) Strain (yp in/in) 
10 12 S4O00 . 192 
30 20 57600 196 
60 2u | 61200 200 
120 28 64800 208. 
180 32 68400 212 
240 32 72000 ~ 220 
300 36 75600 224 
360 © 36 79200 228 
420 36 82800 232 
U8 0 | 36 86400 24.0 
540 40 90000 Zu. 
600 _ 36 93600 | 248 
900 36 97200 292 
1200 | 40 100800 292 
1500 mn 104400 256 
1800 92 108000 260 
2400 96 111600 264 | 
3000 — 60 115200 268 
3600 64 118800 268 
5400 72 , 122400 264 
7200 | 84 : 126000 268 
10800 108 129600 272 
14400 124 | 133200 272 
18000 — 132 136800 276 
21600 144 : 140400 276 
25200 152 144000 _ 280 
28800 156 147600 280 
32400 164 151200 280 
36000 172 154800 288 
39600 y+ eee 158400 292 
43200 180 165600 © 292 
46800 18% 169200 296 


50400 188 172800 te 300 


TABLE A-3 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


Pie Oper. STRESS = 56,000 psi. Compressive 
Time (sec) Strain (p in/in) Time (sec) Strain (p in/in) 

10 32 54000 . 224 
30 U8 57600 228 
60 56 . 61200 232 
120 64 64800 236. 
180 | 72 68400 24.0 
240 68 72000 244 
300 68 75600 248 
360 68 79200 252 
420 co. 82800 | 260 
480 v2 86400 264 
540 76 90000 2712 
600 . £6 93600 276 
900 80 97200 a 
1200 84 100800 | 280 

1500 84 104400 | 288 
1800 88 108000 292 
2400 32 111600 292 
3000 96 125200 296 
3600 36 118800 300 
S400 1o4 122400 300 
7900 > #22 126000 304 
10800 124 129600 308 
14400 LHO 133200 312 
18000 148 136800 312 
21600 156 140400 312 
25200 172 144000 312 
28800 180 147600 312 
32400 188 151200 320 
36000 200 154800 324 
39600 204 158400 324 
43200 208 165600 — 328 
46800 212 169200 332 
50400 220 ~ 172800 7 336 


TABLE A-4 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


ES Boek. STRESS = 66,000 psi. Compressive 
Time (sec) Strain (yp in/in) Time (sec) Strain (yp in/in) 

10 32 54000 . 512 
30 52 57600 520 
60 68 | 61200 528 

120 88 64800 — 532. 
180 96 68400 536 
240 104 72000 548 
300 Li 75600 556 
360 124 79200 564 
420 132 82800 568. 
Ugo | 140 86400 580 
540 144 90000 588 
600 - 148 93600 588 
900 168 97200 600 
1200 | 184 100800 | 608 

1500 196 104400 624 | 

1800 208 108000 632 
2400 224 111600 640 
3000 236 115200 644 
3600 248 118800 652 
5400_ 276 122400 652 
7200 j 296 126000 660 
10800 332 129600 664 
“14400 360 133200. 668 
18000 376 136800 672 
21600 400 140400 676 
25200 412 144000 680 
28800 «428 147600 | 680 
32400 — HUY 151200 692 
36000 4u60 — 154800 696 
39600 472 158400 700 
43200 480 | 165600 712 
46800 488 169200 716 
50400 504 —— 172800 "716 


TABLE A-5 CREEP TEST DATA, ALUMINUM 7075 .T651 (BARE) 


i= “B5°r. - STRESS = 43,000 psi. Tension 
nnn 
Time (sec) Strain (p in/in) Time (sec) Strain (p in/in) 
Peres age Bae eh a cn a eet ee Bae 

60 g 10800 52 

120 15 14400 | 57 

180 16 18000 62 

240 | 18 21600 63 

300 18 25200 68 

360 19 81000 110 

420 Sg 8244.0 123 

480 - 20 82800 124 

540 21 86400 123 

600 20 30000 12 

900 yale | 93600 124 

1200 21 97200 127 

1500 20 100800 ) 130 

1800 ZA. LO¥4O0 130 

2100 23 108000 paral 

2400 24 111600 145 

2700 26 113400 144 

3000 28 163800 142 

3300 29 165600 138 

3600 31 ~ 169200 132 
7200 39 


TABLE A-6 CREEP TEST DATA, ALUMINUM 7075 .T651 (BARE) 


Ser as ee Se. 
Time (sec) Strain (yp in/in) Time (sec) Strain (yp in/in) 
Ci ais A pA TEAS RNS nD Fee ny ce Eee ie 
60 32 25200 110 

120 38 28800 114 
180 LO 36000 120 
24.0 43 39600 12? 
300 | 4 43200 124 
360 | 46 58500 131 
420 48 61200 132 
48 0 : 50 64800 130 
540 51 68400 129 
600 52 72000 129 
1200 oy, 75600 aes 
1800 63 | 79200 128 
2400 65 82800 126 
3000 68 86400 125 
3600 69 93600 125 
7200 ‘78 100800 127 
10800 , 87 108000 130 
14400 93 115200 135 
18000 98 118800. 139 
21600 103 175500 142 


TABLE A-7 CREEP TEST DATA, ALUMINUM 7075.T651 (BARE) 


T ‘eB? F - STRESS = 66,000 psi. Tension 
Ov 
Time (sec) Strain (yp in/in) Time (sec) Strain (py in/in) 
a eae as ee a a eg er Ne eae 
60 42 39600 - d7'S 

600 70 43200 L180 
900 89 46800 185 
L200 68 90400 196 
1800 . 74 94000 202 
2400 78 72000 215 
3000 83 75600 217 
3600 . 90 79200 218 
4200 98 82800 216 
9400 Ls 86400 218 
7200 106 93600 219 
9000 Ja: * 100800 : 224 
10800 | * 223 108000 | ' 22/1 
L4¥4O0 LoS 115200 231 
18000 136 122400 : 236 
21600 141 129600 242 
29200 148 : 162000 256 
28800 Los 165600 209 
32400 160 169200 204 
36000 166 ~ 172800 254 


TABLE A-8 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 
Wy Be Fe STRESS = 72,000 psi. Tension 


mn eooooeemms= 


Time (sec) Strain (p in/in) Time (sec) Strain (yp in/in) 


ot a a a ee SS SSS 


60 12k 61200 826 
120 151 64800 838 
180 173 68400 B54 
240 189 72000 B64 
300 198 75600 872 
360 209 82800 886 
420 27 86400 891 
48 0 224 90000 900 
540 230 93600 907 
600 237 97200 916 
900 260 104400 927 

1200 282 109800 94.3 
1500 297 113400 949 
1800 312 115200 966 
24.00 337 118800 982 
3000 360 122400 998 
3600 377 126000 1010 
5400 u18 129600 1027 
7200 u51 133200 1048 
12180 516 - 136800 1066 
14400 540 140400 1077 
16200 560 144900 1092 
19800 574 147600 1098 
21600 588 151200 1112 
25200 603 154800 | 1120 
28800 627 158400 1131 
46800 756 162000 L140 
50400 773 169200 1140 
54000 790 172800 1145 


TABLE A-9 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


Tet 35°F. STRESS = 43,000 psi. Compressive 
Time (sec) Strain (yp in/in) Time (sec) Strain (yp in/in) 

10 20 54000 372 
30 28 57600 376 
6.0 32 | 61200 388 
120 40 64800 396. 
180 ran 68400 LOY 
2u0 Lg 72000 L16 
300 4g 75600 u20 
360 52 79200 Tae 
420 56 82800 436 
480 60 86400 rarer 
540 60 90000 4u52- 
600 64 93600 460 
900 68 97200 468 
1200 | 76 100800 | 472 

1500 7.6 104400 ugo 
1800 84 108000 488 
2400 92 111600 496 
3000 100 115200 500 
3600 108 118800 504 
5400 132 122400 512 
7200 148 126000 520 
10800 176 129600 524 
14400 204 133200 532 
18000 228 136800 540 
21600 24.8 140400 552 
25200 264 144000 556 
28800 280 147600 560 
32400 296 151200 568 
36000 308 154800 572 
39600 324 158400 580 
43200 336 165600 — 596 
46800 352 169200 604 
50400 356 172800 604 


TABLE A-10 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


10 

30 

60 
120 
180 
240 
300 
360 
420 
480 
940 
600 
900 
1200 
1500 
1800 
2400 
3000 
3600 
9400 
7200 
10800 
14400 
18000 
21600 
295200 
28800 


T 


Time (sec) 


Slo oe es - STRESS 


Strain (py in/in) 


32 
48 
64 
84 
104 
108 
FL 2 
i aie) 
120 
120 
eG. 
a I 
120 
128 
144 
156 
184 
LOZ 
208 
236 
260 
280 
308 
332 
360 
380 
400 


Time (sec) 


32400 
~ 36000 
39600 
43200 
46800 
50400 
54000 
57600 
61200 
64800 
68400 
72000 
75600 
79200 
82800 
86400 
90000 
93600 
97200 
100800 
104400 
108000 
111600 
115200 
118800 
122400 
126000 


= 56,000 psi. Compressive 


Strain (p in/in) 


416 
432 
rarer 
460 
472 
484 
496 
508 
520 
528 
536 
552 
560 
568 
576 
584 
592 
600 
608 
612 
620 
628 
636 
644 
652 
656 
664 


TABLE A-ll1 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


TS  Lsort. ' STRESS = 66,000 psi. Compressive 
nnn ORSON 
Time (sec) Strain (y in/in) . Time (sec) Strain (y in/in) 
eae ee ee an ee ae Rs 
10 40 1500 236 
30 56 1800 252 
60 64 2400 284 
120 - 84 3000 316 
180 104 : 3600 7 340 
240 104 5400 LOY 
300 104 | 7200 456 
360 128 10800 524 
440 | 136 | 14400 576 
480 144 18000 596 
540 148 ° 21600 624 
600 156 25200 . 632 
900 184 28800 644 
1200 224 32400 | 652 


TABLE A-12 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


Ta= absoor STRESS = 30,000 psi. Tension 
en 
Time (sec) Strain (yp in/in) Time (sec) Strain (p in/in) 
a re aa ate eee Fe ce Se a Se 

60. 2 5400 45 
120 6 7200 «50 
180 10 10500 55 
240 13 12600 59 
300 14 14400 62 
360 15 15000 64 
420 17 21600 «68 
480 18 25200 71 
540 19 27000 73 
600 20 82800 104 
900 24 86400 101 

1200 27 90000 100 
1500 28 93600 | 9g 
1800 30 97200 : 93 
2100 31 100800 93 
2400 33 104400 | 85 
2700 34 108000 81 
3000 35 111600 83 
3300 36 115200 81 
3600 38 ~ 169200 106 

172800 100 


TABLE A-13 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


T=: 235°R STRESS = 43,000 psi. Tension 
Time (sec) Strain (p in/in) Time (sec) Strain (p in/in) 
60 12 5400 79 
120 21 6180 8 3 
180 26 6600 83 
240 29 7800 86 
300 31 | 10800 93 
360 34 14400 97 
420 35 18000 101 
480 37 21600 105 
540 | 39 25200 108 
600 39 29400 111 
900 43. 7 84600 161 
1260 50 86400 160 
1500 53 93600 j 160 
1800 57 100800 162 
2400 62 108000 170 
3000 66 115200 IZ 
3600 Ti 171000 188 
4200 75 172800 187 
4800 77 


TABLE A-14 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 
T = 135°F ' STRESS = 56,000 psi. Tension © 


eee eee eeeeeeeeeeeeoaaaounanaonononnnnnmmmmmteeee SSS aos 
Time (sec) Strain (p in/in) Time (sec) Strain (p in/in) 
a haa gas ne er a as ee ag ee ee ee 


60 23 3000 ~107 


120 34 3600 113 
180 40 5400 129 
240 4 8 7200 135 
300 54 | 10800 143 

3360 53 14400 159 
420 53 18000 169 
48 0 55 21600 | i 
54.0 56 77400 236 
600 59 79200 236 
900 70 86400 238 

1200 78 93600 236 

1500 84 100800 241 

1800 90 108000 243 

24.00 101 172800 279 


A-18 


TABLE A-15 CREEP TEST DATA, ALUMINUM 7075 T651 (BARE) 


T = 135°F STRESS = 66,000 psi. Tension 
Time (sec) Strain (yp in/in) Time (sec) Strain (yp in/in) 
60 37 39600 508 
120 68 54900 569 
180 75 57600 574 
2u0 89 | 61200 583 
300 99 64.800 579 
360 106 68400 605 
u20 114 72000 614 
48 0 115) -2 75600 621 
540 124 79200 630 
600 129 82800 638 
900 153 86400 660 
1200 167 80000 669 
1800 193 97200 - 683 
2400 213 104400 703 
3000 231 108000 715 
3600 9h3 115200 » 998 
5400 277 122400 74.6 
7200 301 151200 798 
10800 sk 2 154800 809 
14400 374 158400 805 
25200 U4 O 162000 809 
28800 WEY 165600 814 
32400 477 169200 817 
36000 493 172800 822 


A-19 


TABLE A-16 QUANTITIES ae €, and r FOR ALUMINUM 
7075 T651 (BARE) at 85°F 

stress ° 

(ksi. ) Em et r 
30 Tension (no data) |(no data) (no data) 
30 Compression 267 x 107'| 36.390 | 2.08 x 107" 
43 Tension 591 x 107'| 75.599 | 2.00 x 107" 
43 Compression 214 x 107°] 127.912 | 1.32 x 107* 
56 Tension 100 x 107°| 111.057 42 x it0°" 
56 Compression 635 20 | ek | Be eo 
66 Tension 455 x 107] 147.821 | 2.35 x 107" 
66 Compression 936. x o>? 396.292. ~05 xX i 
72 Tension 300 x 10°°| 601.845 .03 x 107" 
72 Compression (no data) | (no data) (no data) 
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APPENDIX B 


LISTING OF SUBROUTINES 
MATCON AND SHIFT 


In using CREEP-PLAST, the user must modify two sub- 
routines to reflect the material being used in the analysis. 
Subroutine SHIFT contains the description of the material 


creep characteristics according to the relation 


ad 7 
= - + 
EG e, 1 e ) ae 6 


as discussed in Appendix A. The listing of subroutine 
SHIFT shown in Figure B-1 reflects the equations as pre- 
sented in Appendix A. 

Subroutine MATCON (Figure B-2) contains a description 
of the stress-strain curve including the tangent modu- 
lus as a function of the effective stress. The stress- 
strain curve used is taken from MIL-HDBK-5 and reproduced 
in Figure B-3. 

For this analysis the stress-strain curve correspond- 
ing to LT tension has been used. Using a quadratic function 
in terms of the effective stress yields the following: 


for :o.< 61,300 psi © = (9.788 ie i 
for 61,300 <-o < 76,000 psi 


LZ 


Jo" 


e = (9.788 x 107°)o + (5.873 * 107 


Og oS 


+ (3.539 x 10 ~-)o 


for o >76,000 psi e = (5.0 x 107 )o 


For the tangent modulus 
for o < 61,300 E. = 10.22 = 10° psi 


for 61,300 < o < 76,000 psi 


ie 


1/(9.788 x 10°) + (1.175 x 10°--)Jo 


doe 


Ee 


+ (1.062 x 10 ~~)o 


and for o > 76,000 psi. 
E = 2.0 x 16° psi 


These equations have been incorporated into subroutine MATCON. 
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As a result of a coding error in subroutine SHIFT (not 
shown) a slightly different equation for the tangent modulus 


as a function of stress was used. In particular 


o < 76,000 psi. EE, = 10.22 x 10° psi. 
o > 76,000 psi. EY = 2.00 x 10° psi. 

This is a bilinear curve slightly "higher" than the LT- 
tension curve of Figure B-3. The difference between this 
bilinear curve used and the LT-tension curve is negligible 
(< 5%) as far as the effective stress-effective strain coor- 
dinates at the start of continuous compressive load is concerned. 
The formulation of the LT-tension curve is retained in this 
report in order to show clearly the technique'used to represent 


the tangent modulus as a function of the effective stress. 
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APPENDIX C 


FINITE ELEMENT MODEL OF ALUMINUM 
CREEP “TEST SPECIMEN 


The finite element grid used in the CREEP-PLAST analysis 
is shown in Figure C-1 and C-2. This model is one-quarter of 
the aluminum specimen since symmetry has been used to decrease 
the model size. Boundary conditions are discussed in the text. 
Detailed input data corresponding to the grid is included in 
Appendix D. Included in this section is an additional listing 
of the nodal point coordinates for the convenience of the 


user. 
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APPENDIX D 


LISTING OF DATA INPUT FOR 
FINITE ELEMENT CODE CREEP=-PLAST 


This appendix presents the input data used for executing 
the CREEP-PLAST code. Notice that nodal data is printed ina 
format that does not include all the significant digits of the 
coordinates. However, Appendix C contains: the card images for 
the nodal points and elements from which exact nodal coordi- 
nates can be found. The output presented here is for the 
initial run using CREEP-PLAST. 
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APPENDIX E 


PARTIAL LISTING OF COMPUTER OUTPUT 
FROM CREEP-PLAST: 


In this appendix is presented the output from the code 


CREEP-PLAST corresponding to Figures. 4 and 5 in Section II. 
Since quite a voluminous amount of output is generated only 
the more important time steps have been included. fhe entire | 
output is available upon request. Note that times in this 
output are "seconds" and not "hours". The output for the 


following time steps have been included in this Appendix: 


Time Step No. | Time (sec.) 
1 | Peo 
32 | 30.0 
43 1,000.0 
ran 2,000.0 
W5 3,000.0 
“6 4,000.0 
L7 5,000.0 
‘48 6,000.0 
4g 7,000.0 
50 8,000.0 WW 
51 9000.0. <2 aw 
52 10,000.0 
62 | 20,000.0 
72 30,000.0. 
64 50,000.0 
105 } 80,000.0 
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